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Abstract 

We study in this paper the Heisenberg antiferromagnet with nearest neighbours 
interactions on the Husimi cactus, a system which has locally the same topology as 
the Kagome lattice. This system has a huge classical degeneracy corresponding to 
an extensive number of degrees of freedom. We show that unlike thermal fluctua- 
tions, quantum fluctuations lift partially this degeneracy and favour a discrete subset 
of classical ground states. In order to clarify the origin of these effects, we have set up a 
general semi-classical analysis of the order from disorder phenomenon and clearly iden- 
tified the differences between classical and quantum fluctuations. This semi-classical 
approach also enables us to classify various situations where a selection mechanism 
still occurs. Moreover, once a discrete set of ground states has been preselected, our 
analysis suggests that tunelling processes within this set should be the dominant effect 
underlying the strange low energy spectrum of Kagome-like lattices. 
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1 Introduction 



Motivated by the search and understanding of new quantum disordered ground states in 
magnetic systems, recent years have seen a renewal of interest on the properties of frustrated 
Heisenberg antiferromagnetic (HAF) systems with a special attention to their possible re- 
lations with superconductivity. More specifically, a considerable amount of work has been 
done on frustrated systems whose classical ground states exhibit an infinite number of lo- 
cal continuous degeneracies. Typical examples of such systems are the HAF models on the 
Kagome lattice in two dimensions (2D) and the HAF models on the pyrochlore lattice in 
three dimensions. The HAF models on the Kagome lattice has been originally inspired by 
the experiments on the 3 He layer on graphite [I], |2[] and also by the structure of the SC r GO 
compound 0. Other compounds whose magnetic ions live on a pyrochlore lattice have also 
been experimentaly studied by neutron scattering [|J . From the theoretical point of view, the 
ground state of the quantum HAF model on the Kagome lattice is believed to be a quantum 
spin liquid. This model has been investigated by exact diagonalizations on finite lattices 
H, series expansions f|, Large- N calculations @] and also with semi-classical approaches 
|8|. [|, [K], [ll]]. More recently, some theoretical works on the HAF model on the pyrochlore 
lattice have pointed out that this model is already disordered at any finite temperature ]12| 



unlike the xlassical Kagome Heisenberg antiferromagnet which chooses coplanar spin con- 
figurations at low temperatures. The low energy spectrum of the HAF model with spin | 
on the Kagome lattice [13] and the pyrochlore lattice [14] have been analyzed by performing 
exact diagonalizations on small clusters. The results are quite fascinating in both cases since 
there is a singlet-triplet gap and nevertheless a large singlet degeneracy which increases with 
the cluster size. It seems then plausible that the huge classical degeneracy should play a 
major role in order to explain these surprising results. 

Systems with many degenerate classical ground-states were often considered as useful toy 
models before studying for instance spin glasses for which both frustration and randomness 



are relevant. It has been then noticed by J. Villain and coworkers [fLq| , and also by E. 
Shender |]I6 |, that switching on any source of fluctuation (typically, either classical thermal 
or quantum zero-point fluctuations) has a tendancy to lift the degeneracy between these 
classical ground-states. Often this mechanism picks ordered ground-states since they exhibit 
more symmetries and therefore encourage larger fluctuations. This is the so-called "order 
from disorder" phenomenon. Over the last ten years, many interesting contributions have 
been dedicated to this topic, in the context of the frustrated square lattice |L7]], the Kagome 



lattice [10, nj, and the pyrochlore lattice ]I2, |14]]. A common trend in all these works is 
the appearance of at least a partial degeneracy lifting within the ground-state manifold, in 
the presence of classical thermal or quantum zero-point fluctuations. In most cases, both 
types of fluctuations favour the same subset of ground-states. This appears to be a little 



1 



surprising, since we shall provide in this paper an explicit example for which they behave 
in a qualitatively different way. But the question of degeneracy lifting is only one part of 
the story, since it is mostly a local analysis in the vicinity of a given ground-state. Maybe 
more crucial to the physical behaviour is the global question of whether the system remains 
close to any of the favoured ground-states or not. This question is closely connected to the 
size of energy barriers separating these pre-selected ground-states. If these barriers are too 
small, either thermal activation or quantum tunneling may delocalize the system in phase- 
space, and thus destroy the selection itself. In fact, evidence now strongly accumulates to 
support such a scenario, either for the Kagome or the pyrochlore lattices, which exhibit 
classical and quantum spin liquid phases with apparently no long-range order. Although 
this global aspect of the selection mechanism is very interesting, there are few systems for 
which a rather complete analytical treatment is available. Therefore, in spite of the fact that 
different systems often exhibit different behaviours, we believe that valuable insight may be 
gained from simpler models which can be studied more easily than more realistic geometries. 

Among these, the Husimi cactus (see Figures 3 and 4) has already provided an explicit 
example of a system with a continuous manifold of classical ground-states where energy bar- 
riers are too low to freeze the system in the vicinity of a given ground-state in the classical 



thermal case, or to prevent dramatic tunneling processes in the T = case []18||. The rel- 



evance of tunneling is clearly demonstrated by comparing the relatively small ground-state 
degeneracy (for spin S = 1/2) which scales as the number of sites and the huge classical 
degeneracy: the number of coplanar spin ground-states is for instance growing exponentially 
with the number of sites. The question we address here is: tunneling between which classical 
ground-states? More precisely: does the quantum zero point energy depend on the classical 
ground-state? We have found that unlike classical thermal fluctuations for which all classical 
ground-states on the Husimi cactus are equivalent, quantum fluctuations do favour the dis- 
crete subset of coplanar classical ground-states. We believe it is one of the clearest examples 
where thermal and quantum fluctuations behave in a qualitatively different way, already 
at the non-interacting spin-wave level. In order to clarify the origin of this effect, we have 
set-up a semi-classical analysis of order from disorder. The main idea is simply to calculate 
the spin-wave spectrum from the classical equations of motion for spins, since the classical 
frequency also determines the quantum zero-point motion energy of a harmonic mode. The 
classical dynamical approach greatly simplifies for instance the separation between fast oscil- 
lations around a classical ground-state and slow drift motions along the classical ground-state 
manifold, which are connected to zero modes. It also provides a very simple understanding 
for the different nature of classical thermal and quantum zero point fluctuations. In the 
context of the Husimi cactus, the equivalence between all the ground-states for thermal fluc- 
tuations is related to the presence of continuous transformations which conserve the total 
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energy and phase-space volume, and therefore the Gibbs measure. An example of such a 
transformation is a rotation of all the spins below a given site on the cactus around the spin 
at this given site. However, conservation of the Gibbs measure is not sufficient to enforce 
conservation of the zero-point motion energy, since hamiltonian dynamics requires also a 
symplectic structure on phase-space (Poisson brackets), which turns out not to be preserved 
by the aforementioned transformations. We shall show that the non-conservation of Poisson 
brackets under these phase-space transformations is due to the fact that they cannot be 
generated (via Poisson brackets) by functions defined on phase-space. This shows that it is 
rather difficult to preserve quantum-mechanically a continuous classical degeneracy which 
is not enforced by a symmetry of the Hamiltonian. We expect that only a discrete residual 
degeneracy survives after introducing quantum zero point motion. 

We believe this analysis contributes to rule out one possible scenario which might have 
been advocated to understand the unusual situation of the spin 1/2 Kagome lattice. The 
generation of very low energy scales in the singlet sector might have been attributed to very 
slow drift motions close to and along the continuous set of classical ground-states. Note that 
this continuous set is not a smooth manifold for either the Kagome or the pyrochlore lattices, 
since for those the number of zero modes depends on the classical ground-state |10| . Although 
simple and appealing, this idea does not work for two reasons. The first one is the large gap 
between the singlet and the spin one sectors, which shows the impossibility to construct long- 
lived wave packets centered mostly around a single classical ground-state. The other reason 
is provided by our analysis: because of quantum zero-point fluctuations, the ground-state 
continuum turns into a discrete set, so that classical drift motions in the vicinity of this 
continuum acquire finite typical frequencies. Therefore, we conclude that a semi-classical 
analysis of the energy spectrum on Kagome-like lattices (of which the Husimi cactus is a 
simple example) should focus on tunneling processes within a discrete set of preselected (via 
quantum fluctuations) classical ground-states. These processes are crucial to induce the gap 
in the spin one sector. The main challenge for the Kagome lattice is then to understand why 
they lift so weakly the degeneracy in the singlet sector. We haven't addressed this crucial 
question here, but we have clarified what happens when we try to quantize a system with 
a smooth manifold of classical ground-states, whose dimension is larger than the number of 
independent symmetries of the Hamiltonian. Either tunneling plays a minor role, but then 
the classical slow variables experience an effective energy landscape which is no longer flat, 
or tunneling is the dominant effect, as we believe is the case for both the Husimi cactus and 
the Kagome lattice. 

Our presentation is organized as follows. Section 2 will present very simple examples 
with two degrees of freedom and a one dimensional degenerate ground-state manifold. This 
will illustrate how order form disorder works in a classical hamiltonian system, provided the 



3 



fast modes are prepared with a finite value of the corresponding action variable. Section 3 
is dedicated to a simple spin chain for special values of coupling parameters which produce 
a more symmetrical ground-state manifold (the full sphere S 2 ) than one would expect from 
the symmetries of the Hamiltonian (rotations around a single axis). Section 4 presents the 
main results for the Husimi cactus, and section 5 contains our geometrical interpretation 
of the difference between thermal classical and quantum zero-point fluctuations. A brief 
conclusion is given in section 6, and several appendices present some technical details that 
most readers will prefer to skip. 



2 Some simple examples with two degrees of freedom 

In this section, we consider simple examples of dynamical systems where the classical ground- 
state manifold is one dimensional. This will demonstrate explicitely that classical dynamics 
exhibits the phenomenon of order from disorder provided the fast modes are "fast enough" . 



2.1 An integrable toy model 

Let us consider a dynamical system with two degrees of freedom, with coordinates (Xi, Pi, X 2 , P 2 ) 
on phase space. Here (Xi,Pi) denote a slow mode which is continuously degenerate along 
the line (X%, 0, 0, 0) in phase space, and (X 2 , P 2 ) a fast mode, which oscillation frequency Q 
is a smooth function of X 1 . We thus introduce the following Hamiltonian: 

H = \Pl + \n{Xi){Pi + Xl) (2.1) 

We note that Pi commutes with H only if Q is a constant. So the classical degeneracy 
is not in general due to a symmetry of H. This Hamiltonian is integrable since we can 
perform the following canonical transformation (Xi, P±, X 2 , P 2 ) — ► (X±, Pi, J 2 , 2 ), where 
X 2 = y/2J 2 sin 8 2 and P 2 = ^2J 2 cos #2- With these new canonical variables, H becomes: 

H = \Pl + n(X{)J 2 (2.2) 

The motion of the fast mode is described by: 

h = (2.3) 

e 2 = fi(^i) (2-4) 

So, as we expect, the oscillation frequency is controlled by the "slow" variable X%, and is 
therefore time-dependent in general. Here the action J 2 is exactly conserved, which is a 
special feature of this model. In general, J 2 may be approximately conserved with good 



accuracy, since it is an adiabatic invariant [O, [ZOfl. However, this requires a good separation 
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between the slow and the fast time scales. This question will be addressed soon. Turning 
now to the "slow" variables, we have: 



X 1 = P 1 



(2.5) 
(2.6) 



We see that if J 2 is non- vanishing, the slow variables (X 1? Pi) acquire a non-trivial dynamics, 
with an effective potential given by: V e ff(X 1 ) = J 2 Q,{X-y). Therefore the classical degeneracy 
is dynamically lifted, and X\ oscillates around the value which minimizes the fast frequency 
Q(Xi). In a semi-classical approach, the fast oscillator is excited with quantized values for 
J 2 = h(n + 1/2), where n is a non- negative integer. So the effective potential corresponding 
to quantum zero-point motion is obtained from J 2 — \. Although very simple, this model 
captures the main idea we need throughout the rest of the paper. 

2.2 A more generic example 

The special status of the previous model comes from the fact that action-angle variables 
corresponding to the fast motion do not depend on the slow variables (Xi, Pi). In quantum 
mechanical words, only the frequency of the fast motion depends on Xi, but the eigenstate 
basis of the corresponding harmonic oscillator doesn't. In such a case, the adiabatic princi- 
ple is exact, and eigenstates of the time-dependent problem are given by those of the time 
independent one (up to phase factors). But in general, we expect the local eigenstate basis to 
depend on the actual value of X\. As a result, the classical action is only approximately con- 
served, and therefore transitions to excited states may occur for the fast oscillator, specially 
if the two frequency scales are not very well separated. We show here a simple example 
where the adiabatic principle may break down. We think this is a very interesting prob- 
lem in classical mechanics, which deserves further investigation. But this would be beyond 
the scope of this paper. We just provide here an illustration that complicated things may 
happen. Let us then consider the following Hamiltonian: 



This represents as before a potential landscape similar to a gutter with variable curvature 
depending on Xi. Assume first that Xi can be frozen. The fast oscillator evolves with a 
frequency Q(Xi) = (1 + tX\ ) x l 2 . The corresponding action-angle variables (J 2 ,0 2 ) would 
now be given by: 



H=\Pl + \Pl + \{l + eXl)Xl 



(2.7) 



^ = (#T)) 1/2 ^2 

P 2 = {2J 2 Q(Xi)f 2 cos6 2 



(2.8) 
(2.9) 
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But now, since the definition of ( J2, 62) involves explicitely X\, the corresponding transforma- 
tion is no longer canonical. As before |P 2 2 + |^ 2 (^i)^l = ^(-^1)^2 is the effective potential 
of the slow variables, but J 2 is not exactly conserved. As usual, this effective Hamiltonian 



is obtained by averaging over one period for the fast variables ||19|| . To evaluate when the 
adiabatic approximation makes sense, let us compare the slow and the fast frequencies. The 
averaged motion of the slow variables is given by: 

X x = P x (2.10) 
P i = _J 2^(Ai) = -^2e (1+ J£ )(1/2) (2.H) 

For the small oscillations of Xi, we find the slow frequency u = (Ji^) 1 ' 2 which has to be much 
smaller than the minimum value of fl(Xi), that is we require J 2^ <C 1. We have performed a 
numerical integration of the equations of motion for this model in order to check the overall 
picture. The results are summarized in Figures la, lb, lc, Id. The Figures la and lb represent 
the "quasi orbits" of the slow and fast oscillators for two different values of the parameter 
e = 0.1 and e = 0.25. We clearly see a qualitative change in the dynamics as e is increased 
with a breakdown of the adiabatic approximation for larger values of e. This can be seen in a 
more explicit way by plotting directly the variations of the adiabatic invariant J2 of the fast 
motion during one quasi period of the slow motion. This has been done in Figures lc and Id. 
The adiabatic approximation is no longer valid for e = 0.25. It would be very interesting to 
identify the underlying mechanism by which this breakdown occurs. We just mention here 
that Figures lc and Id suggest that J 2 acquires a time dependence with frequencies which 
are close to simple harmonics of the "fast" frequency. Therefore, a perturbative approach 
might give valuable insights here. 

The data support however the validity of our picture up to rather large values of J 2 e, 
around 0.1. We shall now move on to spin systems with many degrees of freedom. 



3 A ferromagnetic chain 

3.1 Ground state and dispersion relations in the different phases 

We study in this section the spin dynamics of a one dimensional XXZ chain with an 
anisotropic quadratic term. The Hamiltonian reads 

n = Y,J (s?s? + i + s?s? +1 + ASfSfa) + E(^) 2 (3.1) 

i i 

where Si is a three dimensional vector normalized to S (the length of the vector is a constant 
of motion from a classical and quantum point of view). 
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The classical phase diagram of this model can be determined easily and has been drawn 



on Figure 2 (see also ref. |2T[] for a quantum analysis with S = |, S = 1). It is composed 



of three different phases: the ferromagnetic one, the planar one and the antiferromagnetic 
one (AF). In order to compute the low energy spectrum in the various phases (at large S), 
it is easy to use the Holdstein-Primakoff transformation. Nevertheless, we prefer to work 
here with classical spin dynamics at T = in order to give a consistent treatment of the 
order from disorder phenomenon. This assumes a spin S reacts like a rotor which can be 
parametrized by a vector S. The evolution of a quantity O(S) is given by: 

dO(S) 



dt 

with {• • ■} the Poisson brakets defined as: 

{/,»} = 



{0,H} 



c abc d f gg S c 



dS a dS b 

with e abc the completely antisymmetric tensor. If S is a vector living on a sphere S 2 , the 
corresponding classical Poisson brackets are given by: {S a , S b } = e abc S c . After a standard 
quantization process, we would obtain the usual commutation rules for spin operators (see 
|22|1 for a rigorous treatment a la Dirac). The equations of motion then become: 

1T~ e os~ bb {6 - 2) 

These equations of motion can also be derived directly from a variational principle, the latter 
being more adapted for a direct path integral quantization scheme. 

We are now ready to derive the equations of motion following eq. ( |3.2| ) for the hamiltonian 
( p.l| ) in each of the three possible phases. 

In the ferromagnetic ground-state, we use the decomposition S = S z z + S 1 - (with 
(z) 2 = 1) and we obtain the system of equations: 



dt 



dt 



J 



(*Sfj~-i + S n+1 ) A S^z + A(S^ l _ 1 + S* l+ i)z A S. n 
D 



+ J (S* n (zAS^) + (zAS^S, 



(3.3) 
(3.4) 



In the vicinity of the ferromagnetic ground-state, we can assume: 

S n (t) = Sz + expi(kn - ut)e + 0(e 2 ) 
with e.z = 0. By solving this system, we find the dispersion relation: 

\u(k)\ = 2JS\A + — -cos(A;)| 

'J 



(3.5) 
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We would have obtained the same low energy spectrum (with H = 1) by using for example 
the Holdstein-Primakoff transformation. When A + -y = —1, the mode at k = tc becomes 
slow and has the expected k 2 dispersion for a ferromagnetic ground-state. 

In the planar region, we use the following general decomposition for S n = S® + SS n 

S n = (-l) n Su + (A(-l) n v + Bz)expi(kn - oot) (3.6) 

where S® = (—l) n S. We have also defined u.z = 0, (u) 2 = 1 and v — zAu. The equations 
giving A and B are: 

-icoA = 2JS[1 + Acos(fca) + 
iujB = 2JS(l-cos(ka))A (3.7) 

As |A| < 1 + in the planar region, we find: 

to = ±2JS{1 - cos(ifea))3(l + Acos(A;a) + (3.8) 

Two interesting limits can be taken. First, when A + y — ■> — (1) + , we recover the soft mode 
in k 2 , whereas when A — — — > (1) _ , we find a linear dispersion around k = typical of an 
AF ground-state. 

In the AF area, we can use the decomposition: 

S n = S z n + = (-l) n Sz+ [(-l) n e + ffi expi(kn - cot) (3.9) 

with e.z = ff.z = 0. By inserting this decomposition in the equation of motion, we obtain 
the following dispersion relation: 

co = ±2 J S '[( A -^) 2 - cos 2 (Jfea)]* (3.10) 

In the limit A - § -> 1+ we recover the same dispersion relation as A — -j — > 1 found in 
the planar region. 

3.2 Order by disorder along the transition lines 

In the following, we take special interest in the transition lines between these areas. Let 
us first analyze the transition between the ferromagnetic and planar ground-state. Along 
the line separating both regions in the phase diagram, the classical ground-state is a priori 
strongly degenerate and can be parametrized by the two spherical angles (6, 0) (6 is the 
angle between S and z). The Hamiltonian (|3.1| ) is invariant by rotation around the z axis 
and therefore the energy does not depend on cf). We are here in the unusual situation 



where the classical ground-state manifold is invariant under a larger group (SO (3)) than 
the Hamiltonian itself (SO (2)). This suggests that quantum fluctuations may reduce this 
larger group to the smaller one by selecting one orbit on the classical ground-state manifold 
under the smaller group. In our view, this is then a toy model for the order from disorder 
phenomenon. We use the decomposition: 

The equations of motion reduce to: 



da 



n 



dt 
dt 



-J 



-J 



K-i + ft + i) A + I A| + a* n+1 ) A ft\ (3.11) 
{(ati + ft + i) A ft] (3.12) 



where we have used A+ -y + 1 = 0. We are looking for solutions around c?o = S(cos8z+sm6u) 
with u.u = 1 and u.z = 0. We therefore write a = a + 5a where 5a can be parametrized 
by: 

5a = [Av + P(cos 9u — sin 9z)\ exp i(kn — ujt) 
with v = z Au. After standard manipulations we find: 

oj e (k) = ±2JS[cos(ka) - 1][1 + (|A| - l)sin 2 #]5 (3.13) 

with a dispersion in k 2 for the slow mode at ka = 0. We see that when |A| ^ 1 there is a 
selection mechanism. When | A| > 1, the spin wave zero-point motion favours a ferromagnetic 
ground-state (9 = 0) whereas for |A| < 1 the planar ground-state is favoured (9 = |). 

This is the most simple example of order from disorder phenomenon in so far as the 
variable 9 has no dynamics (S z = S cos 9 is a constant of motion) and plays rather the role 
of a parameter than a true variable. If we try to classify various possible situations, we may 
say here that we have a dynamical system with one conserved quantity and a 2-dimensional 
ground-state manifold. Furthermore, in the present model, this ground-state manifold is 
generated by two symplectic flows which do not commute (for instance the rotations along 
z and the rotations along x). A general model for this is given by the Hamiltonian: 

N N 

H = \ E E <k3(Pi)PiPj + bijiPjXiXj + 2c ij (P 1 )P i X j (3.14) 

i>2 j>2 

For this Hamiltonian, P\ is conserved and therefore plays the role of cos#. The ground-state 
manifold is given by Pj = Xj = for 2 < % < N. This 2D manifold appears as a collection 
of ID orbits under the symmetry generated by Pi, and characterized precisely by the value 
of Pi. Note that the drift motion along these orbits arises from the excitations of the "fast" 
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variables (2 < i < N), as may be seen from the equation giving the time-derivative of X\. 
Apart from this drift motion, we can clearly separate the fast variables from the collective 
ones (i — 1). 

Returning to our example, it is straightforward to describe the average drift motion 
< cf) > due to the quantum zero-point motion of the internal (i.e. k ^ 0) modes. The idea 
is, as explained before, to view the zero-point motion energy of the internal modes (which 
depends on 9) as an effective potential for the collective variables and 9. Summing up the 
zero-point motion energy obtained from |3.13| over all k's yields the effective potential: 



I , rf (0) = 2JSJ(1 + jsin 2 fl) (3.15) 



This effective potential for S z leads to an average drift for the <fi variable: 

• Dcos9 



1 + 7 sin 2 9 



(3.16) 



Let us now consider the Ising-like case |A| > 1, for which 9 = is selected The most 
interesting fact is that </> still has a non zero limit when 9 — > 0. Indeed 

Um6 = D (3.17) 

How shall we interprete this fact ? It simply means that the rotor has a residual energy 
even when 9 — > which corresponds to an easy-axis ferromagnetism. This is nothing but an 
anisotropy gap of the system. It is easy to check using Holstein-Primakoff 's representation 
that the lowest excitation corresponds to a magnon with a gap D for ka = ±7r. This 
simple calculation gives a classical image for a gap which has however a purely quantum 
origin, through zero-point motion of internal modes. The results ( |3.13| ) ( |3~T7j ) show that 



classical hamiltonian dynamics provides in this case a simple way to understand the low-lying 
excitations of the quantum system. This analysis suggests that the critical theory controlling 
the behavior of the ferromagnetic chain can be described by a simple rotor (independently 
of the nature of the spin S). This has been analyzed in detail and from a different point of 



view m 



We can perfom the same work for the transition line separating the planar ground-state 
and the AF one. The details of the algebra have been gathered in appendix A. The low 
energy spectrum is given by equation ( |A.14| ): 




u e (k) = ±2JS(ka)\l + — sm 2 e (3.18) 



with y = A — 1. For A = 1, the spectrum does not depend on 9. When A < 1, we find 
that spin waves favour the planar ground-state whereas for A > 1, the AF ground-state is 
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selected. We have therefore the same kind of order from disorder phenomenon as for the 
Planar-Ferromagnetic transition. 



Through this example, we have analyzed a simple one dimensional example of order from 
disorder. Similar situations may occur in 2D, frustrated systems like for example in the 



J\ — J2 XY or Heisenberg model on a square lattice [17|, ^4|, where the classical degeneracy 
is global, and partially enforced by a true symmetry. This list is obviously not exhaustive and 
many other models with frustrating interactions experience this kind of order from disorder. 
However, note that the spin chain analyzed here is special (and actually simpler) in that 
the classical ground-state manifold is not isotropic for the symplectic form which defines the 
classical dynamics (see section 5.1) 

The Hamiltonian (|3.1| ) has been studied only from the point of view of order from disorder 
using a semi-classical analysis. A full quantum analysis demands a more careful study for 
instance of tunneling processes which are sensitive to the actual value of the spin (integer or 
half- integer) . This is not the purpose of this article. We refer the reader to [21| for a careful 
study of the phase diagram and transitions. 



4 The Husimi tree 

In this section, we will see another non-trivial example of the order from disorder phe- 
nomenon, which we shall study in an extensive way. We consider a pseudo-lattice called the 
Husimi cactus. A three-generation cactus has been represented on Figure 3. It consists of a 
succession of triangles connected only by their vertices. The geometrical dual of the Husimi 
cactus is a Cayley tree with constant coordination number. The interesting fact is that at 
least locally this pseudo-lattice looks like the Kagome lattice. However, by contrast to the 
Kagome lattice, there is no closed loop of connected triangles. Each node of the cactus 
has one Heisenberg spin. We consider only nearest neighbour AF interactions, namely the 
Hamiltonian reads: 

<i,j> 

with Si =1. An essential feature of the Husimi cactus is the huge degeneracy of its classical 
ground-states. This is because minimizing the energy on each triangular plaquette yields 
the constraint: 

E & = (4.2) 

ieA 

which does not define a unique ground-state configuration. Indeeed, suppose one spin on 
a given plaquette is fixed (say Si), then its two neighbours have to live on a half-cone 
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centered on —S\ with an opening angle of |. The ground-state manifold of a finite cactus is 
then a smooth manifold (i.e. the number of zero modes is independent of the ground-state 
configuration). To describe this manifold, let us assume that the top spin is fixed. Then, we 
have to choose a rotation angle for the two "offspring" spins around this ancestor. Because 
of the hierarchical structure of this lattice, a similar freedom exists at each successive level of 
the hierarchy, until we reach the bottom boundary triangles which we shall call the leaves. Of 
course, we have a complete freedom to choose the value of the top spin, which adds two more 
global zero modes. For a p-generation cactus with 2 P — 1 sites, the phase space manifold for 
the Heisenberg model is a product of TP — 1 spheres S 2 , whereas the ground-state manifold 
is isomorphic to the product of S 2 (choice of the top spin value) by a torus of dimension 
(2(f -1 ) — 1) (i.e. one angle per spin which is not on the bottom boundary). Hence the ratio 

Number of zero modes 2^ p_1 ) + 11 . 

p»l. (4.3) 



Number of degrees of freedom 2(2 p — 1) 4' 

The Husimi cactus with Heisenberg spins is interesting so far as it exhibits a macroscopic 
number of zero modes and we expect spin waves to be very soft. Hence, it is not a priori 
obvious whether order from disorder can occur in such a system. In particular, are copla- 
nar ground-states favoured as in the Heisenberg Kagome antiferromagnets? We emphasize 
however that the degeneracy lifting mechanism is slightly different for "real" systems such 
as the Kagome or the pyrochlore lattices since for them, the most relevant effect amounts to 



maximizing the number of zero modes [12 



4.1 Thermal fluctuations 

As we explained in the Introduction, we should distinguish the question of degeneracy lifting 
within the ground-state manifold with the more global question of actual selection or freezing 
in the vicinity of a particular ground-state. Regarding the first question, it is easy to show 
that classical thermal fluctuations do not lift any degeneracy. The reason has been given 
in the Introduction: a rotation of all the spins below a given site on the cactus around this 
ancestor spin preserves the phase space volume element (since the Jacobian matrix is then 
triangular and rotations preserve volumes), and also the total energy. Therefore such trans- 
formations preserve the Gibbs measure. Clearly any two ground-state configurations may 
be connected by a finite number of such transformations. Our statement then follows easily. 



Although very simple, this point has not been noticed in ref. |Tj|. This paper emphasized 
the other aspect, namely that energy barriers are too small to allow freezing (for instance 
after a quench from high temperatures) in any of the classical ground-state configurations. 
Consequently, we do not expect an order from disorder phenomenon induced by thermal 
fluctuations in this model. This system is too soft for a selection mechanism to work. Other 
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situations with extensive entropy and inequivalent ground-states, where there is a partial de- 
generacy lifting within the ground-state manifold but the global selection mechanism fails, 



have been presented in ref. [FS|. More recently, the study of the Heisenberg model on a 
pyrochlore lattice led to the same conclusions [Hj . Here, maximizing the number of zero 
modes favours colinear magnetic states, but the effective barriers are too weak for the system 
to remain confined in the vicinity of these states. 
How similar is the quantum case? 

4.2 Quantum fluctuations 
4.2.1 General considerations 

As in the previous sections, we wish to understand the classical spin dynamics at T = 
on the Husimi cactus. Again, the main feature of this sytem is its huge degeneracy and 
the fact that these ground-states build up a smooth manifold. If we examine the vicinity of 
the ground-state manifold, we see essentially two kinds of motion, those associated to finite 
frequency oscillations around the chosen ground-state and those associated to what we may 
call slow drift motions. Typically, this can be described by an approximate Hamiltonian 



(defined around the chosen ground-state) which takes the general form |25 



1 1 n 

H = -{Pl + --- + Pf) + - E + 9?) (4-4) 



2 y 1 r/ 2 . . . , 



The (Pi, Qi), 1 < i < r + s are associated to the slow drift motions and (Pj, Qj), 

r + s + l<i<n describe the "fast" oscillators. Note that the (Pi, Qi), r + l<i<r + s 

are associated to global symmetries. The equations of motion are: 



Qj — Pj 
Qj = 
Qj = u jPj 



Pj = 0, l<j<r 

Pj = 0, r + 1 < j < r + s (4.5) 
Pj = —UjQj, r + s + l<j<n 



The degenerate manifold is generated by the vectors associated to Q\, ■ ■ ■ , Q r+S ; P r +i, • ■ ■ , P r +s- 

Our main goal is to determine if the spectrum depends on the ground- state we consider 
and more precisely if zero-point fluctations lift at least partially the classical degeneracy. In 
order to compute the oscillator spectrum, we will proceed as follow: 

(i) First we have to find in phase space the submanifold associated to the drift motions. 
They are defined by initial conditions producing a speed parallel to the ground-state sub- 
manifold. For the general Hamiltonian (|4.4j ), the drift space is generated by the directions 
associated to (Pi, Qi) 1 < i < r + s. 
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(ii) Second, the oscillator submanifold is obtained by taking the symplectic orthogonal 
of the drift motions submanifold. 

By symplectic form, we mean the antisymmetric bilinear form which is defined by G(X, X') 
Y,i=i(PiQi - QiPD, where X = (P h Qi) and X' = (P(, QJ) are any two vectors in the 2n- 
dimensional classical phase-space. The main interest of this symplectic form is that the 
Poisson bracket of two functions / and g defined on phase-space is given by: 

{f,g} = Q(Vf,v g ). 

This bilinear form is non-degenerate (i.e. it is impossible to find any non-vanishing vector 
which is orthogonal to any other vector), so for any linear subspace D, the set D 1 of all the 
vectors orthogonal to any vector of D has a dimension dimD 1 - = 2n — dimD. We shall call 
D ± the symplectic orthogonal of D. These notions will be quite useful for the analysis given 
in section 5. But it turns out that even for basic computations, they provide substantial 
simplification as we shall show shortly. 

In the Husimi cactus case, phase-space is simply (S 2 ) N , where S 2 refers to the sphere 
and N is the total number of spins. If we consider only one spin H and two tangent vectors 
dn\ and dn 2 to n, the symplectic form Q is defined by 

Q(dHi, dn 2 ) = —H.{dH\ A dn 2 ) (4.6) 

It simply corresponds to the area of the parallelogram (dHi, dn 2 ). When we have more than 
one site, this definition extends in a straighforward way 

Q(dn, 5H) = — '^2n i .(dn i A dni) (4.7) 

i 

where (dn, 5n) are tangent vectors in phase space. The equations of motion for this spin 
system are: 

5= E (E^)A* (4.8) 

Here A a denotes any triangle on the lattice. With the present form, it is easy to derive the 
linearized equations of motion in the vicinity of a given ground-state configuration rij. The 
spin at site i is Hi + 5Hi, where Hi.SHi = 0. For an equilibrium state, we have J2jeA a — 
for any triangle A a . Therefore, the linearized equations read: 

^r= E (E sn^AfU (4.9) 

a,ieA a jeA a 

Our goal is now to extract the finite frequency spectrum from these equations. Rather than 
remaining abstract, let us apply this strategy to a 3-generation cactus. 
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4.2.2 Spectrum of the 3-generation cactus 

We consider a 3-generation cactus whose sites have been numbered from 1 to 7 (see Figure 
3). At each site, we define a Heisenberg spin rT,. Let us derive the most general drift motion 
around an equilibrium configuration n il<i<7 . We have then ni + n 2 + n% = n 2 + + n 5 = 
n3 + tTq + n 7 = 0. This motion is generated by a global rotation around a vector r, a global 
rotation of vector eh[, a partial rotation of and n 5 around n 2 (of angle 77), and finally a 
partial rotation of and li? around (of angle (). We are looking for tangent vectors Srii 
(ni.Sni = 0) producing a motion parallel to the ground-state manifold. Therefore, we have 
to solve the system (the arrows on vectors will be omitted in order to lighten expressions): 

hi = r A rii = (5rii + 5n 2 + Sn 3 ) A rii 

h 2 = {r + erii) A n 2 = (Srii + 5n 2 + Sn 3 + 5n 2 + <5n 4 + 5n 5 ) A n 2 
hs = (r + erii) A n 3 = (5rii + 5n 2 + 5ris + 6113 + 5uq + 5n 7 ) A n 3 

= (r + erii + r\n 2 ) A n 4 = (5n 2 + Sn^ + 5n 5 ) A n 4 (4-10) 
h 5 = (r + erii + r\n 2 ) A n 5 = (5n 2 + 5n 4 + 5n 5 ) A n 5 
h 6 = (r + en x + (n 3 ) A n 6 = (5n 3 + 5n 6 + 5n 7 ) A n 6 
«7 = (r + e^i + Cn 3 ) A n 7 = (5n 3 + 5n e + 5n 7 ) A n 7 

The fact that we are in the vicinity of an equilibrium state imposes that n§ and are 
linearly independent as well as and n 7 . This fact is used in the resolution of the system 
( |4.10| ) and we thus find that the submanifold associated to the drift motion is locally defined 
by: 

Sni + 5n 2 + 5n 3 = 

Sn 2 + 5n A + 5n 5 = m\ + r\n 2 (4-H) 
5n 3 + 5n 6 + 5n 7 = en x + (n 3 

We call D\ this subspace. The case e = 77 = ( = corresponds to initial conditions with a 
vanishing speed, namely to the tangent space to the degenerate submanifold, which we shall 
denote by D 2 . We have now to take the symplectic orthogonal of this submanifold D\. More 
accurately, we are looking for vectors Ui G (-Di) -1 such that 



i) Vz, Ui.rii = 

ii) V Srii G Di, nj.(5ni A Uj 







(4.12) 



The procedure to construct explicitly this submanifold is not straighforward and is explained 
in Appendix B. Nevertheless, in order to determine the spectrum, we may work with the 
larger submanifold (-D2)" 1 " which has been defined in appendix B. To give some illustration 
of what these subspaces mean, let us return to the "flat" canonical phase-space used in the 



context of the general Hamiltonian 4.4. The drift space D 1 is as we have said generated by 
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the vectors associated to (Pi, Qi) for 1 < i < r + s, and D± corresponds to the finite frequency 
modes (Pi, Qi), r + s + l<i<n. The space D 2 is spanned by the vectors associated to 
Qi, 1 < i < r + s and Pj, r + l<j<r + s. Its orthogonal D 2 contains the direct sum of 
-Df and the subspace of D 2 generated by the vectors associated to Qi, 1 < i < r. As the 
equations of motion [|75] show, D 2 is stable for the Hamiltonian flow. The velocity vanishes 
everywhere on D 2 L f)D 2 , and the image of _D^~ for the linear mapping (Pi,Qi) — > (Pi,Qi) is 
precisely the subspace associated to the finite frequency modes, namely D±. 
Consider now a vector u G D 2 which can be defined by: 

(4.13) 




with + u 2 + iig = 

We now compute the flow associated to u which is naturally defined as 



v 1 = 


(ui + u 2 + u 3 ) Arii = (u\ + uf) A n\ 


= v{ 






v 2 = 


(u\ +U2 + U3 + U2 + U4 + u 5 ) A n 2 = 


(u a 2 +u b 2 + u%) 


A n 2 


= v% + v\ 


v 3 = 


(ui + u 2 + u 3 + u 3 + U6 + u 7 ) An 3 = 


(u% + u c 3 + u\) 


An 3 


= V 3 + V 3 


V 4 = 


(w 4 + -u 5 + u 2 ) A n 4 = u 2 A n 4 = v 4 








^5 = 


(w 4 + m 5 + u 2 ) A n 5 = «2 A n 5 = 








V6 = 


(w 3 + m 6 + W7) A tiq = u 3 A n 6 = v% 








V 7 = 


(u 3 + uq + U7) A n 7 = u 3 A n 7 = Vj 









(4.14) 



where the velocities (t>j = hi) are given by the linearized equations of motion p~9| . By defining 

v 2 = «2 A ra 2! f 4 = ^2 A n 4 ; t>£ = u% A n 5 (4.15) 
v 3 = M3 A n 3 ; t>g = M3 A n 6 ; u 7 = u 3 A n 7 

we check that (D 2 ) x is stable under the Hamiltonian flow. Indeed, we clearly have: 

< + v 2 + v° = ^ + 4 + 4 = v c 3 + v c 6 + v c 7 = 0. 



Moreover, we can prove that v G (Di)- 1 , as we discussed using the Hamiltonian 

To show explicitly this interesting property here, we have just to show that v is orthogonal 
to any tangent vectors Sn satisfying ( |4.11|) . The demonstration proceeds in a pedestrian way. 
The main ingredient involved is that Q(v,8n) involves only the sums of the n field on the 



various triangles, which are precisely given by the right-hand side of equations P"T. A similar 
reasoning is detailed at the end of Appendix B. 

We now compute the frequencies of the oscillators. We therefore consider the square of 
the application giving the Hamiltonian flow: 

flow flow 2 / a 1 r\ 

U > V > —OJ U (4.16) 
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The application closes because of the stability of D 2 under the linearized flow. We obtain: 

(4.17) 

(4.18) 

(4.19) 

We notice that the equations associated to (u 2 , uf) and to (u 2 ,u 3 ) are closed. We have 
therefore to solve two sytems 



US Li-^ 


- (q) b _|_ r,f 


l\ 1 L\ — I Li 2 


A II., 1 q A n „ 

f\ ii>2 \ "3 '\ "'3y 




—U1 2 U 2 


= ( v b + v c- 
\ 2 1 3/ 


A no = (u?; 


A Tio + A 77.9 


A n 2 


—0J 2 U 3 


= ( v b + v f 
\ 2 1 3/ 


A n-i = (u?; 

o \ 2 


A 71<> + !i? A 77.9 


A 77,q 


9 h 
-UJ 2 U b 2 


= v f A n 2 - 


= ((«2 + «i) 


A \ ft 

A n 2 J A 72 2 






= v 2 A n 4 - 


= ((«* + «§) 


A 72 2 ) A 77 4 




-u 2 u\ 


— v 2 A ns = 


= (K + n c 3 ) 


A ra 2 ) A 72 5 




-tjj 2 u\ 


= f 3 A n 3 = 


= ((«S + «i) 


A n 3 ) A 77,3 




-tjJ 2 U% 


= f 3 A n 6 = 


= ((«£ + «!) 


A 72 3 ) A 77, 6 




—U 2 Uj 


= f 3 A n 7 = 


= ((«§ + «§) 


A 77,3) A 72 7 






(4.20) 



~UJ 2 U 2 = —U 2 + (7/3.712)773 + TjUg 

-u; 2 7i3 = -u% + (u2.n 3 )ra 2 + \u a 2 

-u 2 u h 2 = -{u\ + u c 3 ) + ((u* + ^).n 2 ).n 2 ^ 2 

-^ 2 M ^ = + u c 3 ) + {{u\ + uQ. n 3 ).n 3 

In order to solve the system (A), we introduce the vector z = ^\Anl\\ anc ^ wrr te 

u" = pjZ + A n,, 1 < 7 < 3, 

and 

Pi + P2 + P3 = 

. We find for system (A) that the possible values for uo 2 are |, |, 2. We can solve the system 
(B) in the same coordinate system and we find the same eigenvalues. 

We can thus conclude that the eigenvalues of the massive modes for the 3-generation 
cactus are u 2 = |, u 2 = |, uj 2 = 2 and are doubly degenerate. 

Moreover, the spectrum does not depend on the classical ground-state and the order 
from disorder phenomenon does not occur for this system. 

4.2.3 Order by disorder in the p-generation cactus 

We would like to see now if these results also apply to any finite cactus with p generations. 
The general structure of the calculations is as outlined in the previous paragraph. We also 
work within the linear space D 2 whose definition is the natural generalisation of equations 
4.13| . Since D 2 is stable under the linearized flow, we define a linear mapping of D 2 which 
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sends any tangent vector u\ of into the corresponding velocity Vi. An explicit calculation 
of the secular determinant for the p = 4 cactus is presented in Appendix C. The main 
conclusion is that for p > 4, the zero-point energy does depend on the classical ground-state 
configuration. An illustration of this is given on Figure 5 where the zero-point energy is 
plotted along a special curve on the classical ground-state manifold. This curve corresponds 
to the ground-states which originate from a coplanar configuration and are then obtained 
by rotating all the spins below site 2 (i.e sites 2,4,5,8,9,10,11) by an angle 6 e [0,7r] around 
fi2- The zero-point energy is minimal for 6 = irk, (k integer) and maximal for 9 = (k + 
2)71". Therefore, zero-point fluctuations favour coplanar states, and induce effective barriers 
between such states. We remark that all the planar ground-states are equivalent even after 
introducing zero-point fluctuations. This is because the Husimi cactus has a huge space- 
group which preserves both the Hamiltonian and the symplectic structure. For instance, 
consider one triangle a, f3, 7 in the bulk of the cactus (a being the top site of the triangle). 
If we exchange the sites (3 and 7 and also all their descendants, it is clear the lattice is 
preserved. It is also clear that any coplanar ground-state may be mapped into any other by 
a finite sequence of such transformations. Now, Figure 5 also suggests that these degenerate 
coplanar states are connected by rather low effective energy barriers. Indeed, the height 
of the barrier is about 1% of the low energy excitation which is very low. We would have 
thought intuitively to a much higher barrier since we rotate all spin below sites 2. It indicates 
that the effective low energy landscape is rather flat and that tunelling processes will play a 
very important role. On intuitive grounds, we expect these effective barriers to remain finite 
as the system size goes to infinity, since such a rotation around a given ancestor spin induces 
only local corrections to the small oscillation matrix (see appendix E3). 

Another feature of the cactus is that this degeneracy lifting mechanism occurs only in 
the bulk of the cactus but not on the leaves. The 3-generation cactus appears in fact as 
a particular case because it is built only of leaves (the leaves are defined as the bottom 
boundary triangles). We have demonstrated these results for the 4-generation tree, but the 
conclusion still holds for a p-generation cactus (p > 4) by a straightforward induction proof. 
The special case of leaves is addressed in section 5.3. 

This set of results suggests that physical properties of this system will be dominated at 
T = by dramatic tunneling processes within the discrete set of coplanar classical ground- 



states. This will be the subject of further work [26|. 

A simpler example in the same class is what is referred to in the litterature as the 
Delta chain. It consists of a chain of triangles linked by their vertices as shown in Figure 
6. This system has been much studied and the specific heat was shown to have a double 
peak structure as in the Kagome antiferromagnet f27j. Therefore, it was believed that an 
understanding of its properties can shed light on the Kagome antiferromagnet. Theoretical 



18 



interests on this chain have also been enhanced by its experimental realization p3[. Many 



studies have been devoted to its extreme quantum version namely for S = 1/2 spins [27, ^9] 



The ground-state is exactly known to be a dimer state [EH]]. The low lying excitations 



were found to be kink and antikink-type domain walls in the dimer singlet ground-state [^9 

Our semi-classical analyses applies directly to an open Delta chain, and we find that 
coplanar ground-states minimize the spin-wave zero-point energy here as well. 



5 A geometrical analysis of the spin dynamics on the 
Husimi cactus 



5.1 Properties of the flow associated to rotations around a given 
spin 

Let us consider one spin at a site labelled i. The infinitesimal rotation of angle e around 
this spin defines a vector field on phase-space. More accurately, if the phase space is (S 2 ) N 
(with N the number of spins), a tangent vector rhi to the configuration rlj verifies rhi-fli = 
(i G [1, • • • , N]). If j is a descendant of i then rhj = eniAHj else fhj = 0. We will characterize 
and analyze in the following the properties of these vector fields. 

• The first and main property of these vector fields is that the energy is conserved along 
their flow in phase-space. 

• Secondly, we may wonder if these vector fields can be associated to the symplectic 
flow of a fonction defined in phase-space (that would therefore commute with the 
Hamiltonian!). According to the previous section, we know that the answer to this issue 
is negative because a particular order is selected (indeed, if the hypothesis were true, 
it would imply that actual symmetry protects these internal rotations and therefore 
prevents any selection mechanism from acting). Nevertheless, we now prove it in a 
direct way. In a general manner, suppose we have a vector field U l on a manifold and 
gij denotes the symplectic form in any coordinate system. We assume that: 

3F such that IP = g ij djF (5.1) 

It implies that Vz, aj = gijW = d{F, where aj is a one-form. We consider the simple 
case of two spins n\ and f?2 and take the rotation of n2 around n\ of angle e. This 
defines a vector field mi = 0; m2 = erfi A n^. If dn\ and dn<i are respectively tangent 
vectors to n[ and n,2 then: 

ai.drii + ai.dn^ = — ni.(dfii A mi) — fi2-(dn2 A 7712) 

= —en2-(dn2 A (fi\ A n 2 )) = fii.dn2 (5.2) 
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This property comes from 

ai.dn 1 = drFgijU 3 = Q(dHi,dH 2 ) 



where {?(.,.) has been defined in (|4.6f) . From (|5.2|) , it appears that aidn 1 is not an 
exact differential because it is not closed. Indeed, the equivalent of the rotational is 
defined by 

diaj - djcn = bij, (5.3) 
If dni and 5Hi (i G [1, 2]) are two tangent vectors to the point Hi, then 

bijdn l 5n J = dHxSH^ — bH\dH 2 ^ (5-4) 

We can therefore conclude that those flows associated to internal rotations around a 
given spin are not symplectic flows of fonctions defined on the phase-space despite they 
do not change the energy. 

Another interesting property that may be checked concerns the commutativity of the 
flows associated to two different rotations. Suppose, we first perform a rotation around 
a spin % and then a rotation around a spin j which is a descendant of spin i. This is 
not a priori obvious that the two operations commute. However, we have shown in the 
appendix D that this is the case. 

Finally, we want to mention another property of vector fields. Let us consider two 
internal rotations and their two associated vector fields say U and V. Suppose U 
corresponds to a rotation of angle e around one spin Hi and V to a rotation of angle 
r) around one spin Hj. If i and j are not comparable in the lattice hierarchy then 
G(U, V) = 0. But if j is a descendant of i, then Q(U, V) has no reason to be equal to 
0. 

Nevertheless, on the ground-state manifold, it is easy to verify that 

Vf W S| Ground-State^ = ( 5 ' 5 ) 

This follows directly from the fact the ground-state property YjS[ a \ i indicating the 

i 

three sites of a given triangle a. In more technical words, the submanifold of the ground- 
state manifold generated by these internal rotations is isotropic for the symplectic form. 



5.2 Vector flow associated to internal rotations and Hamiltonian 
Flow 

In this section, we want to clarify why the flow associated with these internal rotations does 
not in general commute with the Hamiltonian flow (except for the rotations of the leaves of 
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the Husimi tree). As we shall see, this is closely related to the non- vanishing of the 2— form 



6y 0- 



5.2.1 General considerations 



In this subsection, we derive some general considerations which will be useful for our prob- 
lem. We consider a phase-space with a symplectic closed form Q (of components g^), a 
Hamiltonian H which is invariant under the flow associated to a vector field X. It reads as: 

C X H = X%H = 0, (5.6) 

with Cx the Lie derivative in the X direction. We also define Y to be the Hamiltonian 
vector field which therefore satisfies 

Y l = g i3 djH (5.7) 
The quantity we are interested in is the Lie bracket: 

[X, Yf = X i d j Y i - Y j d j X i (5.8) 



By using the definitions (|5.71 , |5.8|) and the invariance of H under X ( |5.6|) , we can easily prove 
that: 

[X, Y]* = - {-(d j9 lk )X> + g ij (djX k ) + gt^X*)} (d k H) 

= -a ik d k H (5.9) 

where we have defined a tk = —(djg lk )X j + g lj (djX k ) + g jk (djX l ). 
The relation (6.91) can also be formulated as: 



[X,Y] i = -^a jk Y k (5.10) 

with: 

Oij = 9ika M 9ij = -aji = {d k 9ij)X k + gi k (djX k ) + g kj (diX k ) (5.11) 
We can therefore conclude that: 

[X, Y] = <=^> aij Y j = (5.12) 

The 2— form plays a crucial role. We first restrict to the simplest hypothesis Vi, j = 0. 
This condition has a simple interpretation in differential geometry and means that X is a 
Killing vector field ||32|| . It means that the symplectic form ^ is invariant under the flow 
associated to X, i.e. 
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c x g = o 



or in components: 

(d k9ij )X k + g lk {8 3 X k ) + g kj (diX k ) = (5.13) 

Therefore = •<=>- CxG = 0. Moreover, if X is a Killing vector field, it is at least locally 
the flow associated to a conserved quantity F |32] (a conserved quantity because the flow of 
X leaves H invariant |. 

From these generalities, we can conclude that a sufficient condition for the flow U asso- 
ciated to internal rotations to commute with the Hamiltonian flow is that U should be the 
flow of a conserved quantity. We have proved in section 5.1, that it cannot be true for our 
symplectic form. 



Nevertheless, we must not forget that the condition ( 5.12 ) is weaker. In the Husimi 
cactus, we have Oy ^ as already suggested by eq. ( |5.4|) and as will be seen further. A 
natural issue occurs, namely can we find a vector field X leaving the Hamiltonian invariant 
such that dij 7^ but a^Y^ = ? Before analysing this question, let us make a pause to 
summarize the important consequences of the equivalences derived above. 

The relation (|5.12p implies the following statement: 

In a system with a huge continuous classical degeneracy, the order from disorder phe- 
nomenon (induced by quantum fluctuations) will in almost cases occur unless these degen- 
eracies are protected by true symmetries of the Hamiltonian. Note that here, order from 
disorder refers to degeneracy lifting within the ground-state manifold and not to the global 
behaviour of the system. 



5.2.2 Application to the Husimi cactus 

We now turn back to the question where ^ 0. Can we find in the Husimi cactus vector 
fields such that a^Y^ = ? Suppose we consider the vector field X associated to the rotation 
of angle e around of the vectors fij (fij are descendant spins of fii). The 2— form is 
also given by di(gj k X k ) — dj(gikX k ) which corresponds to the 2— form bij defined in (|5.3|) . 
Therefore, if Uj and Vj refer to two tangent vectors to the point ty, then according to relation 

(0, 

a(u,v) = ^^Ui.Vj — Uj.Vi (5-14) 

where the notation j > % means j descendant of i. We now consider the 1-form fa = a^V^ , 
where V° is the Hamiltonian flow at the point fij in phase-space. For the Husimi tree, it is 

§ To show that, we suppose X 1 = g lj djF <==> d L F = gijXK By using the Schwartz equality didjF = 
djdiF and the Leibnitz rule, we prove the equivalence. Note that this also requires that gij is closed i. e. 
di9jk + dj9ki + d k gij = 0. 
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(5.15) 



(3(u) = a(u, v) = ^2< ^(n k A tij) > Ui 

j>i [keNj J 



fceiVi j>i 



Z! {n k /\n i ).Y J U j 



(5.16) 



This summation has absolutely no reason to vanish. We can thus conclude according to ( |5.12 ) 
that the Hamiltonian flow does not commute with the flow associated to internal rotations 
(except the case of the leaves of the Husimi tree as we will see in the next subsection). To 
analyze this phenomenon more precisely, we shall also consider the vicinity of an equilibrium 
state, and study the evolution of the small oscillation matrix under the action of the flow 
X associated to internal rotations. This has been developed extensively in appendix E. We 
only summarize the main results. First, we have recalled in appendix El how the variation 
of the oscillation matrix A under the flow X is connected to the 2— form a«: 



It is another way of showing that the low-energy spectrum is invariant with the flow of X 
when X is a Killing vector field. 

We have then applied this general formulation to the Husimi cactus in appendix E2 and 
proved explicitly that the oscillation matrix is not invariant under the flow associated to 
internal rotations (see eq: ([E.15| )). 

And finally, in appendix E3, we have compared the oscillation matrices for two configu- 
rations connected by an internal rotation and shed light on the general dependence of the 
low energy spectrum with the angle of the rotations. 

5.3 Case of the leaves of the Husimi tree 

We now treat the problem of the leaves of the tree. We consider one leaf (n ,ni,n 2 ) with 
(ni,n 2 ) the two edge spins. We have seen in appendix B that the degeneracy still holds. 
According to our geometrical discussion in subsection 5.1, this degeneracy must therefore be 
associated to a conserved quantity. 

For one triangle, the Hamiltonian simply reduces to 




(5.17) 
(5.18) 




H = (n + n x + n 2 ) 2 



(5.19) 
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This degeneracy still occurs because the quantity (Hi + n 2 ) 2 commutes with the Hamiltonian 
(under Poisson brackets). Indeed, the symplectic flow associated to (Hi + n 2 ) 2 is given by 



^ = n 2 An\ 



§ = Hi A n 2 (5.20) 



§ = if {1,2} 
where A is the parameter associated to the transormation. Therefore, 

d J^tM = o and d 4 = (5.21) 
dX dX 

It is worth noticing that such a transformation cannot be in general identified with all the 

rotations we may imagine, namely: 

-a rotation around n 

-a rotation around n[j excluding H 

-a rotation around H® including n*o 

Here, n° denotes an equilibrium value, and we assume Hi is close to n°. The only exception 
corresponds to (Hi + H 2 ) colinear to no for the first case and to n[] for the last two cases. 

Furthermore, we can prove explicitly that the order from disorder phenomenon does not 
apply for the leaves of the cactus. We will simply check in the following that the flow 
associated to the function / = (Hi + n 2 ) 2 commutes with the Hamiltonian flow. If we replace 
in Figure 7, the application R by /, we have therefore to show (using the same notations) 
that f = f. 

We consider an infinitesimal transformation associated to the fonction / = (Hi + n 2 ) 2 . 
Therefore, 

(Hi, n 2 ) — ► (Hi + en 2 A Hi; H 2 + eHi A H 2 ) 

Let us consider a tangent vector (fji,fj 2 ) to (Hi,H 2 ) in phase-space. The symplectic transfor- 
mation / implies that: 

(Hi+ffi, n 2 +rf 2 ) — > (n 1 +rfi + en 2 An 1 + en 2 Arfi + erf 2 An 1 ; n 2 +ff 2 + eniAn 2 + en 1 Aff 2 + erfiAn 2 ) 

(5.22) 

Therefore 

(rfi, m) — > A; ifa) = (m + ™ 2 Arf 1 + eff 2 A Hi, n 2 + eH x Aff 2 + effi A H 2 ) (5.23) 

— * 

The Hamiltonian flow in the initial configuration defines the ^ i e {0, 1, 2} as follows 

Co = -n§ A + *h + V2 + Vo + ff-i + ff- 2 ] 

g = -n° A [fjo +ffi + f) 2 ] (5.24) 
k £ 2 = -n° 2 A [rfo + ffi + ff 2 ] 



24 



where the sites (—1) and (—2) are the NN of site (0), in the ancestor triangle. 
Because of ip = ff and ijjx + ^2 = Vi + ^2, we have 

fjT = [fjo +ffl + 772] =*pT = [ipO + fa] 



Therefore, following the notations of Figure 7, 

' Xo = -n§ A [$r + % + #-1 + ^-2] 
< xi = -[n? + en" A n$] A r/ T 

X2 = -[^2 + A ^2] A $T 
By applying the inverse transformation of / to the vectors x%i we can easily prove that 



Therefore, the existence of a true symetry (associated to (Hi + n 2 ) 2 ) leads to the absence of 
order from disorder for all the spins belonging to the edge of the Husimi cactus. Note that 
it is in complete agreement with the geometric vision of order from disorder developed along 
this paper. 

6 Summary, Conclusions and perspectives 

Our main conclusions have been in fact already presented in the Introduction. The main 
result of this work is to show that continuous degeneracies which are not enforced by con- 
tinuous symmetries of the Hamiltonian are more fragile quantum-mechanically than under 
the action of classical thermal fluctuations. This is because classical thermal fluctuations 
are not sensitive to the underlying Poisson bracket structure which is of course crucial for 
quantization. To put this work in perspestive, we give in a Table [l] a summary of what we 
believe is the present status of three closely related systems: the Kagome, the pyrochlore 
lattices and the Husimi cactus. The first line of the Table |I| refers to the local analysis 
of degeneracy lifting. We have included both thermal and quantum fluctuations. Note 
that the degeneracy lifting mechanism is quite different between the "real" lattices and the 
Husimi cactus. For the former, optimal states are the ones which maximize the number of 
zero modes. This generates a discrete subset of classical ground-states (namely the copla- 
nar states for the Kagome lattice and the colinear states for the pyrochlore lattice), within 
which a finer selection may occur (for instance favouring an ordered state on the Kagome 
lattice). Note that the zero mode counting mechanism does not introduce much difference 
between classical thermal fluctuations and quantum zero-point fluctuations. However, the 
residual finer selection is sensitive to whether the system is classical or quantum mechanical. 




(5.26) 
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This mechanism is of course not available on the Husimi cactus, for which the number of 
zero modes is independent of the ground-state configuration. The second line of the table 
summarizes the global selection effects namely whether the system remains or not close to 
any of the favoured ground-states. 

Two extensions of the present work may be considered. The first one would be to set-up 
a similar semi-classical analysis for systems in which the number of zero modes actually 
depends on the ground-state configuration. As we have just mentioned, this would be more 
relevant to "real" systems such as the Kagome or the pyrochlore lattices. The other direction 
would be to take advantage of the simpler phase-space structure on the cactus to analyze 
tunneling processes within the discrete set of coplanar ground-states. As usual, these effects 



are sensitive for instance to the value of the quantum spin S ||33|| , and a rich pattern of possible 
quantum phases as 5* varies may occur in such a system. Numerical diagonalizations would 
therefore be very instructive, specially if performed for several values of the spin S. 



Acknowledgements: We would like to thank P. Azaria and F. Mila for several interesting 
discussions concerning this work. 
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Classical Husimi cactus 


Quantum Husimi cactus 


Local degeneracy 


No 


Yes (coplanar) 


lifting 


this work 


this work 


Global 


No 


No 


selection 


[US] 


[18| 





Classical Kagome 


Quantum kagome 


Local degeneracy 
lifting 


Coplanar 

or ordered 10 


Ordered 
1-0 


Global 
selection 


Coplanar states 



No (Spin Liquid) 






Classical Pyrochlore 


Quantum Pyrochlore 


Local degeneracy 


Colinear state 


Colinear state 


lifting 


[34 




Global 


No 


No (Spin Liquid) 


selection 


[12] 


§1 



Table 1: Status of the ground-states of the classical and quantum Heisenberg model on the 
Husimi, Kagome and pyrochlore lattices with respect to local or global selection phenomena 
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A Low energy spectrum along the Planar- AF transi- 
tion 



We suppose -j — A + 1 = with A > 0. Because of the AF behaviour, we use the decompo- 
sition S n = (— l) n a n + al n with a n .l n = 0. a represents the lattice constant. Because of the 
anisotropy in the z direction, we are forced to separate a n and l n in two parts: 

°n = °n + °n I l n = l n + ^ 

We can write the equations satisfied by these four vectors (we directly take the continuum 
limit): 

Ha L -> - 8 2 rr z 

-j- = AaJ^ Aa z ) + Da\l L A^ Y )+4a{D + J){l z Aa L ) (A.l) 
^- = 4aJ(/ ± A<? ± ) (A.2) 

Hp fPrr 1 

As in section 2, we are looking for solutions around S° = S(cos6z + sin8u) with u.u = 1 
and u.z = 0. We write S = S° + 5S where 5S = (— l) n 5a + al. Using the above equations 
we obtain the following system for Sa z , Sa- 1 , l z , I . 



ir ->_\_ 

= 4a J S {l 1 A cos 6 z) +4a{D + J) S(t z A u sin 6) (A.5) 

(JjL 

= 4aJS{l ± Acos dz) (A.6) 

Hp- d 2 Sa z d 2 Sa ± 

= -(^D + J)a 2 S sin 6{-^— A u)-Ja 2 S cos 6(^-4-^ A z) (A.7) 
at ox 2 ox 2 

Hp 3 2 ()cf^ 

u = _j a 2 5sin #(^_ A £) (A. 8) 
at ox 2 

We use the following form for Sa and /: 

5a = [Az A u + Bw] exp i{kx — ut) 
I = [CzAu + Ew]expi(kx -ut) (A. 9) 



where w = — sin 6z + cos 9u. Using the basis (z } u,v = zAu), and the system ( |A.5| - |A.8| ) , we 
obtain four equations: 
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iooB = -4a JC (A. 10) 

iuoA = AE[J + D sin 2 0] (A.ll) 

iuE = -J(ka) 2 A (A.12) 

tuC = (J + D sin 2 9) (ka) 2 B (A.13) 



The resolution of this simple system gives the dispersion relation 



uo = ±2JS{ka)Jl + j sin 2 9 (A.14) 



with j = A — 1. 



B Vectors associated to the spectrum for the 3-generation 
cactus 

In this appendix, we are looking for vectors Ui G (Di) such that: 

i) Vi, Ui.rii = 

ii) V fe» G A, ^2ni.(5ni AUi) = (B.l) 

i 

where D\ is defined by ( |4.11| ). We first determine the symplectic orthogonal of the subspace 
D 2 defined by the solutions of the homogeneous equations: 

5rii + 5n 2 + Sn 3 = 

5n 2 + 5n 4 + 5n 5 = (B.2) 
Sn 3 + 5tiq + 5n 7 = 

We have obviously D 2 C D\. Therefore the symplectic orthogonal of D\ is included in the 
one of D 2 . It is easier to work with the submanifold D 2 since D 2 is defined by an intersection 
of subspaces defined by only one condition: 

m = {<W £ K = 0} 

where a indicates the a th triangle and i G {1,2,3}. Hence D 2 = f]D 2 . Therefore the 

a 

symplectic orthogonal of D 2 is given by: 



2! 
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since the symplectic form is non-degenerate. This elementary algebra has interesting conse- 
quences because we have only to work on one triangle. Indeed, the first thing to notice is 
that a vector m tangent to an equilibrium configuration rtj and belonging to (-Df) vanishes 
on all the sites which do not belong to the triangle labelled by a. To show this, we remark 
that any tangent vector 5rii for which 5n\ = Sn^ = Sn^ = belongs to D?;. So Ui has to 
be orthogonal to all such <5rij's. More precisely, this implies J2i^A a n i-(^ n i A Wj) = for any 
set of 5rii defined for i ^ A a (and of course drii.rii = 0). Since Ui.rii = 0, this implies that 
Ui = for any i (jL A a . Therefore, we have only to determine U\ for % G A a . For this, it is 
convenient to consider a unit vector z a normal to the plane containing (n", n^, n^) (note that 
Y^=i n i — 0). The tangent vectors Srii and Ui to restricted to the plaquette A a are then 
conveniently expressed as: 

Ui = Pi z a + ai{z a Anl) (B.3) 
5nf = X iZ a + fii(z a An?) (B.4) 

After an elementary algebra (using n\.n^ = —1/2 for i ^ j), we find that D% restricted to 
A a is the set of tangent vectors such that A x + A 2 + A 3 = and p,\ — p 2 — Hz- This 
can be shown to imply that p\ + p 2 + p% = and o\ = a 2 = 03. So (D 2 l ) ± is the subspace 
of D?; of tangent vectors which vanish everywhere but on the triangle A a . To summarize, 
and changing notation slightly, for each triangle A a containing the sites a(a), /3(a), 7(a), the 
tangent vectors uf in (-Df) -1 are defined by: 

i) Vz, ui.rii = 

U) u* = Mi± a(a), /5(a), 7(a) (B.5) 
Hi) < (o) + ufo a) + u a ^ {a) = 

One vector of D% is written naturally as itj = u i- 

a 

In order to form D^, we have just to impose the orthogonality to the particular solutions 



of the inhomogeneous equations ( f4.11| ). The symplectic form of a vector u in with a 
tangent vector {Sn} equals: 

^[«>N = -EE^« A N ( B -6) 

a i 

According to ( |B.5| ), u\ = if i does not belong to the triangle a, we can therefore focus on 
one particular triangle a. With the same notations as before, the contribution of the triangle 
a to the symplectic form reads as: 

Q[u a , 6n) = -cr(Ai + A 2 + A 3 ) + p x p\ + P2P2 + pzPz (B.7) 
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It is then easy to show that the symplectic form S[u a , 5n] depends only on the sum 8n\ + 

8n 2 + 5n 3 



g[u a ,5n] = {5nl + 5n a 2 + 5n%). 



(B.8) 



The main point is that the sum 8n\ + 5n 2 + 5n 3 is just the quantity involved in the equations 
fld.llf) . The dimension of (-D2)" 1 " is 9. For (Di) ± , we have three more conditions. Its dimen- 
sion is therefore 6 and we have consequently only three independent massive modes. More 
explicitly, for the 3-generation cactus (see Figure 3 for notations), (-D1) 1 " is parametrized by: 



a , a j a 

p a i+p a 2+pi 

P b 2+P\ + P\ 

pl + p% + Pi 



and they are constrained by: 

■2 



Tlx. 









z b A (p b 2 n 2 + p b 4 n 4 + p b 5 n 5 ) - a b z b 



+ 



m. 

n 2 . 
n 3 . 



z c A (p c 3 n 3 + p c 6 n 6 + p c 7 n 7 ) - a c z c 



z b A (p b 2 n 2 + plm + p b 5 n 5 ) - a b z b 



3 
2 

3 
2 

-z c A (pgn 3 + p%n % + p^n 7 ) - a c 2; c 









(B.9) 



(B.10) 



C Vectors associated to the spectrum for the 4-generation 
cactus 

We consider the 4-generation cactus whose sites and triangles have been labelled following 
Figure 4. We follow exactly the same construction as for the 3-generation cactus. We thus 
define vectors u° i G [1, .., 15]; a G {a, ..,<?} such that 



uf.rii = 

uf = if i is not a site belonging to the a th triangle 

th 



(c.i; 



E 



U.: 



0, i is a site of the a thiangle 



The vectors u define the symplectic orthogonal D2 of the degenerate submanifold D 2 . The 
subspace is stable under the linearized hamiltonian flow. This defines vectors vf as: 



{4 + uDArij , j G {1,2,3} 
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3 

J 3 
„9 



K + ^ + ^An, , j G {2,4,5} 
{u% + u{ + u 9 7 ) A rij , j G {3, 6, 7} 
(4) A n,- , j G {4, 7, 9} 
(u b 5 )A nj , j G {5, 10, 11} 
K) A n,- , j G {6, 12, 13} 

(w?) A % > i e ( 7 > 14 ' 15 l 



(C.2) 



As in the previous case, we can show that the vectors v are in fact in the submanifold D^. In 
order to find the eigenvalues of the flow associated to the vectors u, we consider the square 
of the linearized hamiltonian flow application. We have therefore the following system to 
solve: 





\n% 


+ 4 


+ "s) An 2 + (ti3 + «|+ii 


?)An 3 ] An,, j G {1,2,3} 


(C.3) 




\u\ 




A n 2 + 4 


A «4 + M5 A 725 


An,, j G {2,4,5} 


(C.4) 






+< 


An 3 + Mg 


A n 6 + A n 7 


An,, j G {3,6,7} 


(C.5) 




K 


+ 4 


+ uf) A n 4 


A Tij, j G {4, 


3,9} 


(C.6) 






+ 4 


+ u%) A n 5 


A n„ j G {5, 10, 11} 


(C.7) 




Xa 


+ 4 


+ Uj) A n 6 


An,, j G {6,12,13} 


(C.8) 


2 g 




+ 4 


+ wf) A n 7 


A^-, j G {7,14,15} 


(C.9) 



Let us notice a rather surprising fact: this sytem divides into two subsystems of independent 
variables. The first subset contains equations ( |G.4| , |C.5|) associated to triangles (b), (c), the 
second subsystem contains equations (|C.3|) , (|C.5| - |C.9| ). Therefore, these equations connect 
next nearest neighbour triangles to each other. This property is clearly general and does not 
depend upon the number of generations of the Husimi cactus. This is associated to the fact 
that we consider the square of the Hamiltonian flow application. 

We first consider the systems of equations ( |C.4| , |C.5| ) related to triangle (6) and (c). We 
define a system of coordinates where: 



ni 



;i,0,0); n a = (~ ^? 0); n 3 = (-^, ~ 0), 



2 2 ' v 2 2' 

where z a = (0, 0, 1) is a vector orthogonal to the plane containing n\, n 2 , n 3 and is linked to 
triangle (a). In a similar manner, we can define vectors z\ z c associated to triangles (b) and 
(c). The vectors u decompose on this basis as follows: 



u: 



■If: 



p\z b + a"{z" A n\ 



G {2,4,5} with Ep' = ° 



p)z c + a c (z c A n$) , j G {3, 6, 7} with ^ Pj c = 



(CIO) 
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We first consider the equation ( |C4| ) with j = 2: 

u 2 u b 2 - {u\ + ul) + n 2 .{u b 2 .n 2 + u c 3 .n 2 ) + -{u\ + u\) + n A .{u\.n 2 ) + n 5 .(u b 5 .n 2 ) = (C.ll) 

We then use the properties of the vectors Uj, namely u b 2 + u\ + = and w^.r^ = 0- We 
project this equation on the z An 2 axis and use the relations: 



u\.n 2 



u 5 .n 2 



2^ 



The vectors 2 6 and z c can be parametrised respectively with the angles 9b and 9 C which 
represent the angles made with z a . Therefore, we have 



e fe ^sin(^),e 6 isin(^),cos(^) 



/ c v / 3 . , 
e —y - sm( 



-e - sml 
2 v 



cos 



with e b , e c = ±1. After projection on the z An 2 axis, the equation (|C.11| ) becomes: 

1 



(uj 2 -3)a b + a c 



+ Pt 



^ cos 9b cos 9 C + e e c sin sin 9 C 
I 

—e c cos 9b sin # c — e b sin 9 b cos # c 







(C.12) 



Using similar algebra, we project (U.ll) on the z axis and obtain 



>2 



?3 



+ a c 



cos 6*b cos 9 C + e 6 e c sin 9 b sin C 
e c cos 9b sin C — — sin cos C 



(C.13) 



If we consider equation ( C.4T) with j = 4, 5, we obtain equations compatible with equations 
( p.l2| , [CT3|) . We can treat the system ( |C.5| ) in a similar way and we obtain in this case: 



{uo 2 - 3)a c 



a 



+ P\ 



- cos 9b cos 9 C + e b e c sin 9b sin 9 C 



E b cos 9 C sin 9 b — e° sin # c cos 9b 



(C.14) 



2 §><* 



+ (7 C 



cos 9h cos 6> c + e b e c sin 9 b sin t 



e cos 9 C sin # b sin 9 C cos #6 

2 



(C.15) 
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In the coplanar case (9b = 9 C = 0), we obtain the following eigenvalues for u 2 : ~, |, | with 
double degeneracy. When (9b, 9 C ^ 0), we have to compute the eigenvalues of a 4 x 4 matrix. 
The characteristic polynomial reads: 

9 



P(w 2 = X, 9i, 9 2 ) = X 4 - 9X 3 + 28X 2 - (^36 + - (cos 2 9 b + cos 2 9 c )j X (C.16) 

175 45 63 9 

H sin 2 9 b sin 2 # c H (sin 2 9b + sin 2 # c ) cos 9 b cos C sin 9 b sin 6* c 

16 16 16 4 

and therefore depends on the ground-state ! 

Notice that this result is in apparent contradiction with the case of the 3-generation 

cactus where the spectrum does not depend on the ground-state. In the 4-generation cactus, 

a selection mechanism holds. It is easy from (|C.16|) to see that the configuration minimizing 



the energy is the planar one (9b — 9 C — 0). Under these conditions we find four eigenvalues 

2 = 1(1) 5(2) 7W 
W 2 ' 2 ' 2 ' 

Therefore, the spins belonging to the triangles (a), (b), (c) are coplanar. Rather than 
solving the heavy system of equations ( |C.3|) , ( |C.6[ - |C.9| ) , we can use an extension of the result 
proved in the previous appendix. Indeed, if we consider one centered triangle b connected 
to three free triangles by its vertices. The results of appendix B extend in a straighforward 
way, and no selection mechanism is found in this system with 9 spins. Therefore, the three 
pairs of spins can rotate independently. If we apply this result to triangles (b), (a), (d), (e) 
((b) being the centered triangle), we find that the spins (8, 9) and (10, 11) are free to rotate 
though the spins of triangle (a) are fixed. Consequently, to summarize, there is no selection 
mechanism in the leaves of the Husimi tree! 

D Lie brakets between vectors fields associated to in- 
ternal rotations 

In this appendix, we check in a pedestrian way that the flow associated to two different 
internal rotations commute. In the Husimi tree, we consider Ri(9) and Rj(r]), with Ri(9) 
the rotation of angle 9 around ?v Three cases have to be looked upon. We have therefore to 
compare the result of the application Ri(9)Rj(rf) with Rj(rf)Ri(9) on a vector n k (k ^ i,j). 

• Hi and j are not connected by a descendence relation, the result is trivial. 

• If k is between i and j, then 

R i (9)R j (ri)(n k ) = Ri(9)(n k ) 

R j (r ] )R i (9)(n k ) = Ri(9)(n k ), (D.l) 

which ensures the result. 
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• If k is below both i and j, the result is a priori less trivial. Let us assume that i is an 
ancestor of j. We have to compare: 



R[n h 6]R[n h ri)]{n k ) and R[R{n iy 9)(fij);r]]R[ni, 8}{n k ) (D.2) 
where R[n, a] denotes the rotation of angle a around n. 

A little algebra shows that this property indeed holds and that the two quantities 
defined in ( p.2|) are equal. 



E Vicinity of the equilibrium manifold 

E.l General point of view 

We first analyze the vicinity of a sub-manifold £ of equilibrium state of an Hamiltonian H . 
Therefore, £ is defined by: 

f H(x) =0 

xG^ n \\ E.l 

\ m( x ) = o 

Let xq G £. We consider a small deviation around xq and thus write x l = x l + y l with 
<< ||xo||. The equation of motion for y is: 



^■ = g ij (x +y)d j H(x + y) (E.2) 



Using \\y\\ << \\xq\\, this equation becomes: 

g ij (x )d j d k H(x )y k + O(y 2 ) (E.3) 



dt 

Defining A 1 -{xq) = g^{xQ)djd k H{xo) 1 the linearized equation of motion around x° simply 
reads: 

% = A'jWv* (E.4) 

It is worth noticing that A 1 • is not in general a tensor. It can be checked easily that djd^H 
does not transform in a covariant way under coordinate transformations. However, if we 
restrict ourselves to the submanifold £, the matrix A 1 Axq) does define a tensor. Let us now 
express the invariance of A under the flow associated to X on the submanifold £ : 

C X A = (9 fc A i J )X fc + A i h d j X h - A k 3 d k X l = (E.5) 

with £ the Lie derivative. 

Let us now extend this definition in the vicinity of £, though A is no longer a tensor. We 
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therefore define = (d k A i j )X k + A\djX k - A k j d k X i . 
Using A\ = g^djdf-H, we obtain 

M\ = («9 fe /)X fc («9^tf ) + g a X k {d k d l 8 3 H) + g a (d j X k )(d l d k H) - g kl (d^^d^H) (E.6) 

The invariance of H under the flow of X is expressed by CxH = or in coordinates 

X k d k H = (djX k )(d k H) + X k (d j d k H) = (E.7) 

(dfiX^idkH) + (djX^dkH) + (diX^djdkH) + X k (d J d k d l )H = (E.8) 

The term (djdiX k )(d k H) equals zero on £. By injecting ( |E.8|) in ( |E.6| ), we find that: 

M\ = {(d k g u )X k -g ik (d k X l )-g kl (d k X i )}(d l d 1 H) 

= -atdtdkH, (E.9) 

where a u has been introduced in (15.90. The relation (IE. 91) can be still written as 



= -g im a mn g nk g k[ (d l d J H) = -a ik A k j (E.10) 

Not surprisingly, we recover the 2— form corresponding to Cxg- We can therefore con- 
clude that if X is a Killing vector field, then the matrix A of small oscillations around an 
equilibrium configuration is invariant under the flow associated to X. 

E.2 Application to the Husimi cactus 

Let u and v two tangent vectors to an equilibrium configuration {n } of phase-space. Ac- 
cording to (|E.1CQ . 

M(u, v) = -u i a ik A k j v j = -a{u, Av) (E.ll) 
We recall how to build the matrix A for the Husimi tree. We consider Hi = n? + Vi with 



Hi.Vi = 0. Therefore, 



( A fyi = -jr = E n 3 A fk 



^ n°A^+ ^ ^ An° + 0(w 2 ) (E.12) 



where Ni indicates all the nearest neighbours of i. At lowest order, 



(Av), = w { = -n° A [ Yl v k + \ztVi ] (E.13) 



where is the coordination number of site i. Using the expression (|5.14|) of a(u, w), we find 



a(u, Av) = )XujWj — UjWi) (E.14) 

j>i 
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where w has been defined in the previous relation. The expression 

12 ™3 = - 12 H °j A ( 12 Vk + \zjVj 
j>i j>i \keNj 

can be simplified by noticing that: 

• For k > i and k £ N i: the term containing v k vanishes because ( I] n° + \z k f^ = 0) 
for an equilibrium state. 

• For k > i and k G N iy the term containing v k is: 



12 n°j + \z k n\ j Av k =n° i Av k 

\3&N k 1 J 

• For k = i, the term containing Vi is n° A Vi 
From these three types of contributions, we infer: 

a(u,Av) =n°. | ^ v k + v { I A ^ i + ftj. I ^ u fc + 2^ I A i (E.15) 
[ \k>i,keNi J J [ \fceiVi / j>i J 

Consequently, the 2— form is clearly different from zero. This result proves that rotations 
of spins fij around (with j descendant of i) do not preserve the form of the oscillation 
matrix around equilibrium. Notice also that it explains why a direct calculation of the low- 
energy spectrum with standard methods is difficult to implement. However, this result does 
not enable us to show explicitly that the matrices A 1 ■ of two equilibrium configurations are 
not equivalent. In order to prove it, we have to compute the characteristic polynomial and 
show it depends on the chosen ground-state as we checked explicitly in appendix C for the 
4-generation cactus. It seems difficult to get general analytical expressions as in appendix 
C for a p-generation cactus. Another less ambitious solution consists in analyzing how the 
oscillations matrices of two different equilibrium configurations are connected. This is the 
subject of the next subsection of this appendix. 



E.3 Links between two different classical configurations 

Consider one vector r] in phase-space. We want to compare the Hamiltonian flow of rj with 
another one linked to 77 by a rotation around n°. Furthermore, it provides us with a way 
of connecting the matrices of small oscillations around two different classical ground-states. 
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This is schematized in Figure 7. We are to compare £ and £. In this sheme, we first begin 
with a state 77. The state ip is deduced by a rotation around n° of angle 9, therefore, 

i> 3 = R{f%, efa ifj>i 
ipj = rjj elsewhere 

We have to distinguish the cases j > % and j > i. We first examine the Hamiltonian flow 
around this new equilibrium configuration. It defines a state x as follow 

Xj = -i?(n°,0)[n°]A [J2 ^k + \z^A if j > i (E.17) 
Xj = A \J2 + \ z i^3 ] elsewhere (E.18) 



We have used ( E.13|) in order to establish these equations. We now apply the inverse rotation 



R(n^, —9) (see Figure 7). We also use the property discussed in appendix D: 

Rin^^RiR^^e^^Rin^e) = R^rf) (E.19) 
In the limit 77 — > 0, we deduce: 

R{n\- -9) {i?(n°, B)[n!j] A R{ff i7 9)[u}} = ff j Au (E.20) 
We have to distinguish four different cases: 

• j > i and j has not % as a nearest neighbour (NN). In this case 
Therefore, 

h = —n°. A I V m. + 
j > i and j NN of i. In this case ipi = rji. Using ([E.20 ), we show 



& = A ( E % + \wi I = (E.21] 



& = -^ A ( E % + J R(n°;-% + ^ % ] (E.22) 



J = I 



& = -a?A[ E E + ^ ] < E - 23 > 

\keNj =i ,k<i k£Nj =i ,k>i 

Notice that the notation k < i means that k is not a descendant of i. 
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• j ^ % and j < i. In this case, we trivially show that £j = £j. 

We may also consider the case of the rotation of the n/s around n° for j > i, and similar 
conclusions can be drawn. 

Consequently, we can infer from this analysis that the problem of small oscillations around 
a new configuration (deduced by an internal rotation around spin n°) reduces to a problem of 
small oscillation in the initial configuration but with a modification of the bounds between % 
and its two immediate descendants (the new "interaction") depends of 0. This modification 
seems to be the origin of the dependence of the low energy spectrum with the angle 0. Note 
that after this change of basis, the modification of the linearized equation of motion due to a 
change in 9 is only local. This is why we do not expect large effective energy barrier between 
different coplanar ground-states. 
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FIGURES CAPTIONS 



FIG. la: Poincare sections of the slow (1) and the fast (2) motions for e = 0.1 with rescaling 
of the variables (X2, P2). The adiabatic approximation seems good even for rather high 
values of e. 

FIG. lb: Poincare sections for the slow and fast motions for e = 0.25. We begin to see the 
deviations from the circular orbits signaling violation of the adiabatic approximation. 

FIG.lc: Variation of the adiabatic invariant J2 of the fast oscillator during one "period" of 
the slow oscillator for e = 0.1. This curve is more accurate to indicate the deviations of the 
adiabatic approximation. 

FIG. Id: Same as figure lc but with now e = 0.25. Deviations from adiabatic motions are 
stronger. 

FIG. 2: Classical pha se diagram of the anisotropic Heisenberg spin chain corresponding to 
the hamiltonian ( |3.1|) 

FIG. 3: The 3-generation Husimi cactus. 

FIG. 4: The 4-generation Husimi cactus 

FIG. 5: The energy barrier separating two discrete ground-states in the 4-generation cactus. 
Spins below site 2 (i.e. at sites 4, 5, 8, 9, 10, 11) are rotated around by an angle 9, giving 
a smooth path on the ground-state manifold. Note that 6 = and 9 = it correspond to 
coplanar ground-states. The height of this barrier is around 1% of the low energy excitation 
which suggests a quite flat low energy landscape in the spectrum. 

FIG. 6: The so called Delta chain with 5 triangles. 

FIG. 7: Schematic representation of the action of the linearized Hamiltonian flow and the 
flow associated to internal rotations. 
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